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Abstract. - We calculate the W ± 'yj + X-production cross sections at next-to-leading order QCD 
for Tevatron and LHC collisions. We include leptonic decays of the W to light leptons, with all 
off-shell effects taken into account. The corrections are sizable and have significant impact on the 
differential distributions. 
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Introduction. — At hadron colliders such as the 
CERN Large Hadron Collider (LHC) and the Fermi- 
lab Tevatron, electroweak boson production in associa- 
tion with jets represents important signal processes as 
yell as backgrounds to future searches beyond the Stan- 
dard Model (BSM). One example is the measurement of 
anomalous tri-boson couplings, arising from BSM physics, 
which can be obscured by higher-order QCD effects. For 
these searches, significance-improving strategies include 
jet-vetos, which amount to subtraction of a leading order 
cross section [1], and are plagued at present with typical 
QCD-scale uncertainties. Improved QCD-precision of pro- 
duction cross sections is therefore essential and has been 
agreed on as a common goal of precision phenomenology 
in the so-called "Les Houches wish-list" [2] . Considerable 
progress in completing this task has been accomplished, 
cf. [3-7] . Concerning electroweak boson production in as- 
sociation with a jet, the QCD-corrections to W + W~ + jet 
have been recently provided in [5] . 

In this letter we examine l / F ± 7+jet production at next- 
to-leading order (NLO) QCD, including leptonic decays of 
the . We devote special care to the development of a 
fully-flexible, numerically stable parton-level Monte Carlo 
implementation, based on the Vbfnlo- framework [8]. Al- 
though full leptonic decays of the massive are in- 
cluded, we will refer to the processes as W ±r yj production 
in the following. 

Elements of the Calculation and Checks. The 

leading order contribution, at 0(a 3 a s ), to the process 
pp — » l~i>t"ij + X includes subprocesses of the type 
qQ — > £~D£"/g, and qg and Qg initiated subprocesses which 




Fig. 1: Representative Feynman graph contributing to the vir- 
tual corrections to the partonic subprocess ud — » e~u e yg at 
0{a 3 a 2 s ). The crosses mark other points where the photon is 
attached to the quark line and the W boson. 



are related by crossing. 

The 10 Feynman graphs of each subprocess can be clas- 
sified into two categories: First, configurations where the 
photon is emitted from the W or the W's decay lep- 
ton, and, second, graphs where the photon is emitted 
from the quark line. Performing the virtual correction at 
0{a 3 a 2 s ), these topologies give rise to self-energy, triangle, 
box, and pentagon (sub-) diagrams. The loop corrections 
are treated using standard methods: Self-energy, trian- 
gle, box and pentagon integrals are evaluated in terms of 
tensor coefficients [9, 10] in dimensional reduction, after 
having applied MS-renormalization. We combine the vir- 
tual corrections to groups that include all loop diagrams 
derived from a Born level configuration, i.e. all self-energy, 
triangle, box and pentagon corrections to a quark line with 
three attached gauge bosons are combined to a single rou- 
tine. This method leaves us with a universal set of virtual 
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Fig. 2: Sample Feynman graph contributing to the partonic 
real emission subprocess ud — > e~u e fgg at 0(a 3 al). The gluon 
is attached to the quark and gluon lines at positions marked 
by the circles. Feynman graph topologies, where the photon is 
radiated off at different positions analogous to fig. [TJ are not 
shown. 



building blocks, which are then assembled for the specific 
process under consideration. This strategy has already 
been applied to various phenomenological studies at NLO- 
QCD precision, e.g. [6,11]. 

The reduction of the loop diagrams has been calculated 
in two independent ways for verification reasons. The 
first approach uses in-house routines within in the frame- 
work of FeynCalc [12] and FeynArts, while the second 
one relics on FeynArts, FormCalc, and LoopTools 
[13,14], with modifications, in particular to the treatment 
of divergencies, as described in [15]. We find that both 
calculations numerically agree within Fortran precision 
for different points. Performing the NLO-computation 
in the chiral limit, the arising infrared (IR) singulari- 
ties have been determined separately in independent ap- 
proaches, and checked against existing results in the liter- 
ature [16,17]. 

The IR singularities encountered in the real emission 
contributions are regularized using the Catani-Seymour 
dipolc formalism [18]. The numerical implementation of 
the dipoles has been numerically checked against Mad- 
Dipoles [19]. The code is optimized such that intermedi- 
ate dipole-results are stored and reused in order to avoid 
redundant calculations. Remaining finite collincar terms, 
after rcnormalizing the parton distribution functions ac- 
cording to [18], were analytically calculated in two inde- 
pendent ways. We integrate the finite collincar terms over 
the real emission phase space by appropriately mapping 
the LO-phasc space, as done in [20]. The cancellation 
of virtual IR singularities against the one-parton phase 
space-integrated dipoles' has been checked analytically. 

We evaluate the leading-order matrix element, as well 
as the subtraction terms using partly modified Helas- 
routines [21] generated with MadGraph [22]. Due to the 
increase of subprocesses when going to the evaluation of 
the IR-subtracted real emission matrix element, optimiza- 
tion is imperative in order not to jeopardize CPU time. 



Table 1: Comparison of integrated W~^j and W~^jj tree- 
level cross sections at the LHC. The cross sections were cal- 
culated with our modified version of Vbfnlo, MadEvent 
v4.4.21, and Sherpa v.1.1.3. The QCD-IR-safe photon- 
isolation is replaced by a conventional separation Rj~, > 1 for 
all jets. We also require Re~, > 0.4 and Rjj > 0.7. All other 
parameters and cuts are chosen as described in the text. 



Here, the matrix element is calculated using the spinor 
helicity formalism of [23], and intermediate numerical re- 
sults, common to all subprocesses, are stored and reused, 
thus speeding up the numerical code. The real emission 
matrix elements, cf. fig. for sample graphs of the par- 
tonic subprocess ud — > e~D e ~/gg, have been checked nu- 
merically against code generated by MadGraph for ev- 
ery subprocess. Integrated results were checked against 
Sherpa [24]. Tab. [1] representatively gives the result of 
our comparison of integrated cross sections with MadE- 
vent v4.2.21 and Sherpa v.1.1.3 for the leading order 
and the real emission dijet contribution, i.e. the process 
pp — > e~D e jjj + X at 0(a 3 O!g), for cuts specified below. 

Concerning the Monte-Carlo implementation of the vir- 
tual corrections, we have implemented the loop contri- 
butions using our Vbfnlo routines, that involve the 
Passarino-Veltman reduction scheme [9] up to boxes, the 
Dcnncr-Dittmaier reduction scheme [10] for pentagons, 
and the spinor helicity formalism of [23]. Throughout, the 
numerical integration is performed using a modified ver- 
sion of Vegas [25], which is part of the Vbfnlo package, 
with different channels for the two- and three-body decay 
of the W boson. Finite width effects of the W boson are 
taken into account using a modified version of the complex 
mass scheme of [26]: The weak mixing angle is taken to 
be real, while using a Brcit-Wigncr-propagator for the W 
boson. This scheme corresponds to the implementation in 
MadGraph. 

For a more detailed discussion of the calculation and its 
numerical implementation, we refer the reader to a sepa- 
rate paper [27]. 

Numerical Results. — We use CTEQ6M parton dis- 
tributions [28] with a s {mz) = 0.118 at NLO, and the 
CTEQ6L1 set at LO. We choose m z = 91.188 CcV, 
m w = 80.419 GcV and G F = 1.16639 x 10~ 5 GcV" 2 
as electroweak input parameters and derive the electro- 
magnetic coupling a and the weak mixing angle sin 6 W via 
Standard Model-tree level relations. The center-of-mass 
energy is fixed to 14 TeV for LHC and 1.96 TeV for Teva- 
tron collisions, respectively. We only consider decays 
to one family of light leptons, e.g. W~ — > e~v e , and treat 
these leptons as massless. The CKM-matrix is taken to 
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Fig. 3: Comparison of the scale dependence of the total cross 
section of pp — ► e~u e ^fj+X at LO (dashed), NLO-QCD (solid), 
and NLO-QCD with the second jet vetoed (dot-dashed) for the 
cuts chosen as described in the text at the LHC. 



Fig. 4: Comparison of the scale dependence of the total cross 
section of pp — > e + v e r yj + X at LO (dashed), NLO-QCD (solid), 
and NLO-QCD with the second jet vetoed (dot-dashed) for the 
cuts chosen as described in the text at the LHC. 



be diagonal, and we neglect bottom contributions through- 
out. A non-diagonal CKM matrix decreases our leading- 
order LHC-result at the per mill level as gluon-induced 
processes dominate the cross section. The correction for 
the Tevatron results, which are mostly quark-induced, is 
about 3%. These corrections are well below the residual 
scale dependence at NLO-QCD. The bottom contributions 
are negligible and can be further suppressed by b-tagging. 
Jets are recombined via the kx algorithm [29] from mass- 
less partons of pseudorapidities \rj\ < 5 with resolution 
parameter D = 0.7. The jets are required to lie in the 
rapidity range \yj \ < 4.5 with pj? > 50 GeV. The photon 
and the charged lepton are chosen to be rather hard and 
central, p e T > 20 GeV, p£ > 50 GeV, |?y 7 | < 2.5, while 
being separated in the azimuthal angle-pseudorapidity 
plane by R h = (A^ + A^) 1 / 2 > 0.2. For the separa- 
tion of the charged lepton from observable jets, we choose 
Rij > 0.2. A naive isolation criterion for the partons and 
the photon spoils IR-safety, yet isolation is necessary to 
avoid fragmentation contributions. We apply the method 
suggested in [30], demanding 



^ ^ ^parton,2 ^ 



i,Ri^<R 



1 — COS R 
1 — COS Si 







Pt 



VR < S , 



(1) 



where the index i runs over all partons, found in a cone 
around the photon of size R. For the cut-off parame- 
ter, that determines the QCD-IR-safe cone size around 
the photon, we choose 8$ = I. 

At leading order, we find a QCD scale dependence of ap- 
proximately 11% for W ± jj production at the LHC, when 
varying /i^ = by a factor two around 100 GeV, cf. 
figs. and [3] for identified renormalization and factoriza- 
tion scales. This scale dependence is only reduced to about 
7% when including NLO-QCD precision for W ± jj. This 
is due to the renormalization scale dependence of the di- 
jet contribution at NLO. Vetoing additional jets results in 
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Fig. 5: Tevatron-comparison of the scale dependence of the 
total cross section of pp — » e~v e ^j + X or pp — > e + v e r yj + X 
at LO (dashed) and NLO-QCD (solid) for the cuts chosen as 
described in the text. 



a stabilization of the cross section, as the veto projects 
on true VF ± 7J events. This agrees with the results on 
W + W~j production [5] and W ± j production [1,31]. 

The difference of W +r fj compared to W~"fj is predom- 
inantly due to the different parton distribution functions 
of the dominant subprocesses. Qualitatively, the findings 
of the W~jj channel generalize to W + jj, accompanied 
by an overall increase of the cross section of about 54% 
(see also Tab. HJ). 

At the Tevatron, fig [5j we find a LO scale dependence 
of 23%, which is reduced to about 8% at NLO-QCD. 
A jet veto is not necessary to stabilize the perturbative 
corrections as additional jet radiation is sufficiently sup- 
pressed by the hard cut on the jet transverse momentum, 
> 50 GeV. 

For the scale choice fi = 100 GeV the total NLO re- 
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Fig. 6: Differential distribution of the photon-lepton separation 
Re~f at LO (dashed) and at NLO (solid). The lower panel shows 
the differential A-factor. The dotted line denotes the total K- 
factor of Tab. fj 
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Table 2: Next-to-leading order cross sections and iGfactors 
for the processes pp — > W^j + X at the LHC for identified 
renormalization and factorization scales, [in = fj,p = 100 GeV. 
The cuts are chosen as described in the text. 



suit differs by about 41% for W~jj from the total LO 
cross section. As usual, however, the total X-factor, de- 
fined to be K = <7 NLO /<7 L , reflects only partly the impact 
of the QCD-quantum corrections on the entire processes' 
characteristics. Quantitative understanding thereof can 
be gained from differential A'-factors of (IR-safe) observ- 
ablcs O, 



K(O) 



d<7 NLO /da LO 



dO 



dO 



In figs. [S] and [71 we representatively show the lepton- 
photon separation and the pr-spectrum of the hardest 
jet at LO and NLO, accompanied by the respective dif- 
ferential A-factors. The distributions develop signifi- 
cant changes when including NLO-QCD precision, yield- 
ing large relative modifications around the total A factor. 

Summary and Outlook. — We have presented first 
results on the NLO-QCD corrections to pp — > W^jj + X 
and pp — > W 7J + A, including leptonic decays and full 
off-shell effects for the W boson. The calculation has been 
implemented in a parton-level Monte Carlo program based 
on the Vbfnlo framework which, thus, is fully flexible ex- 
cept for the limitation that the Frixione definition of pho- 
ton isolation as given in Eq. ([T]) must be used. Using this 
program, we give sample results for total next-to-leading 
order cross sections, as well as differential distributions 
and differential A-f actors. 

We find a fairly reduced scale dependence of the total 
cross sections, cf. figs. for a fixed scale choice = 
[iR = fi and our cuts. The corrections turn out to be 
sizable, around 41% for W~^j production and 34% for 
W + 7j production at the LHC. The total correction at the 
Tevatron is about 30%. 

These total corrections are accompanied by significant 
modifications of up to 60% for differential distributions 
when going from LO to NLO, figs. [6] and [Jj 

A more detailed investigation, including analysis of 
the impact of anomalous couplings and the calculation 
of NLO-QCD jet veto efficiencies for searches suggested 
in [1], is underway Eventually, this process will be made 
publicly available as part of the Vbfnlo package. 



Fig. 7: Maximum jet-pr distribution at leading order (dashed) 
and next-to- leading order QCD (solid). The dotted line de- 
notes the total if-factor of Tab. 
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